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The  interfaeial  energies  of  twin  boundaries  and  staeking  faults  in  metal  magnesium  have 
been  ealculated  using  first-principles  supercell  approach.  Four  types  of  twin  boundaries  and  two 
types  of  stacking  faults  are  investigated,  namely,  those  due  to  the  (1011)  mirror  reflection,  the 
(1011)  mirror  glide,  the  (1012)  mirror  reflection,  the  (1012)  mirror  glide,  the  II  stacking  fault, 
and  12  stacking  fault.  The  effects  of  supercell  size  on  the  calculated  interfaeial  energies  are 
examined.  The  calculated  interfaeial  energies  are  85.5,  81.0,  118.1,  120.0,  8.1,  and  21.8  mJ/m  , 
respectively,  for  the  six  types  of  atomic  configurations. 
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Magnesium  alloys  are  increasingly  being  used  in  a  wide  range  of  applications  due  to  their 
light  weight  and  high  strength.  One  of  the  current  research  frontiers  on  Mg  alloys  is  to 
understand,  estimate,  and  improve  their  low  plastic  formability  to  operate  under  increasingly 
demanding  conditions.  At  the  atomic  scale,  the  plastic  formability  is  closely  related  to  the  ease  of 
the  formation  of  planar  defects  along  the  close  packed  planes,  namely,  twin  and  stacking  faults 
[1-5].  Experimentally,  direct  measurements  of  interfacial  energies  require  very  delicate 
techniques  and  the  results  often  show  large  uncertainties.  For  instance,  the  measured  stacking 
fault  energies  [6-8]  (and  references  therein)  for  magnesium  are  scattered  in  the  range  of  50  -  280 
mJ/m  .  Alternatively,  the  steady  improvement  in  both  computer  power  and  the  efficiency  of 
computational  methods  in  the  past  few  years  has  enabled  the  calculation  of  defect  energetics  at 
reasonable  computational  cost.  [9]  For  the  case  of  metal  magnesium,  existing  first-principles 
results  are  scattered  among  specific  types  of  twin-boundary,  stacking-fault,  approximation  to 
exchange-correlation  energy,  and  supercell  size.  [10-14]  The  main  purpose  of  this  paper  is 
report  the  calculation  of  interfacial  energies  of  twinning  and  stacking  faults  in  magnesium  from  a 
unified  theoretical  framework.  In  particular,  we  report  our  calculated  results  for  the  (1011) 
mirror  reflection,  the  (1011)  mirror  glide,  the  (1012)  mirror  reflection,  and  the  (1012)  mirror 
glide  twin-boundaries  together  with  the  II  fault  and  12  stacking-faults.  The  effect  of  supercell 
size  on  the  calculated  interfacial  energies  will  be  discussed. 

Interface  supercell:  Since  we  employ  first-principles  approach  with  periodic  boundary 
conditions,  the  interfaces  due  to  twinning  and  stacking  faults  are  modeled  using  a  supercell.  In 
fact,  the  crystallographic  theory  of  twinning  [2]  is  rather  complicated  for  a  hexagonal  close 
packed  (hep)  metal.  For  the  special  cases  of  (1011)  and  (1012)  twins,  following  Morris  et  al. 
[11],  the  supercells  are  built  by  first  transforming  the  hep  structure  into  the  orthogonal  structure 


(see  Fig.  1)  by  whieh  both  the  (1011)  and  (1012)  planes  ean  be  seen  more  clearly.  Secondly,  the 
a  lattice  vector  of  the  supercell  is  taken  as  that  parallel  to  paper  surface  and  within  the  (1011)  or 
(1012)  plane  and  the  b  lattice  vector  of  the  supercell  is  taken  as  that  vertical  to  paper  surface  and 
within  the  (1011)  or  (1012)  plane.  Furthermore,  the  c  lattice  vector  of  the  supercell  is  derived  in 
the  direction  perpendicular  to  both  a  and  b.  Then,  cut  away  the  atoms  that  above  the  interfacial 
plane  that  pass  through  a  and  b,  followed  by  making  a  mirror  reflection  with  respect  to  the 
interfacial  plane  that  pass  through  a  and  b.  The  last  step  is  to  shift  the  atoms  that  are  nearest  to 
the  interfacial  plane  to  the  interfacial  plane  for  the  mirror  reflection  twin  or  to  shift  the  reflected 
atoms  by  b  /2  for  the  mirror  glide  twin.  The  examples  of  the  built  supercells  for  the  (1011)  and 
(1012)  twins  are  illustrated  in  Fig.2  and  Fig.  3,  respectively.  The  supercells  for  the  II  and  12 
stacking  faults  are  relatively  easier  to  build.  Take  the  case  of  32-atom  supercell  (16-layer)  as  an 
example,  they  are  just  the  ABABABABCBCBCBCB  and  ABABABABCACACACB 
arrangements,  respectively,  for  the  II  and  12  stacking  faults.  [12]  Note  that  by  supercell  approach, 
each  supercell  contains  two  interfaces.  The  effect  of  supercell  sizes  on  the  calculated  interfacial 
energies  is  studied.  The  built  supercells  with  different  sizes  are  listed  in  Table  1. 

Computational  settings:  To  calculate  the  0  K  energies,  we  employed  the  Vanderbilt 
ultrasoft  pseudopotential  [15]  within  the  generalized  gradient  approximation  (GGA)  [16]  to  the 
exchange-correlation  energy  as  implemented  in  the  VASP  package.  [17,  18]  The  plane  wave 
energy  cutoff  is  132.7  eV  which  is  an  automatic  value  when  setting  the  key  “Prec=High”  in 
VASP.  The  calculation  of  interfacial  energy  involves  finding  the  difference  between  two  total 
energies.  We  have  therefore  adopted  very  dense  A:  points  (See  Table  1)  and  the  Gamma  centered 
scheme  is  used.  We  extracted  the  interfacial  energy,  by 


C  ~  (.^Supercell  ^/icp)/2-'4, 


noting  that  E super  cell  is  total  energy  of  the  supereell,  Eju-p  is  total  energy  of  hep  Mg  sealed  to  the 
supereell  size,  and  A  is  interfaeial  area  whieh  is  sealar  value  of  the  eross  produet  of  lattiee  veetor 
a  and  b.  We  note  that  in  the  above  proeedure  of  extraeting  the  interfaeial  energy,  the  elastic 
energy  resulted  from  the  formation  of  interface  has  not  been  considered  separately.  This  is  a 
reasonable  approximation  for  the  twin-boundary  in  hep  Mg,  as  it  will  be  seen  that  the  calculated 
interfaeial  energies  are  just  weakly  dependent  on  the  supereell  size  (See  Table  1). 

Results  and  discussions:  The  calculated  interfaeial  energies  are  summarized  in  Table  1, 
together  with  and  the  measured  stacking  fault  energies  [6-8]  and  those  previously  calculated 
within  different  approaches.  [4,  10-14]  For  the  twin-boundaries,  it  is  observed  from  our 
calculations  that  the  effects  of  the  gliding  of  the  interfaeial  crystal  planes  on  the  interfaeial 
energy  are  minor,  and  the  interfaeial  energies  of  the  twin-boundary  are  one  magnitude  larger 
than  those  of  stacking-faults.  It  is  further  seen  the  effects  of  the  supereell  size  on  the  calculated 
interfaeial  energies.  For  the  twin-boundary,  increasing  the  supereell  size  up  to  80  atoms  can  only 
decrease  the  interfaeial  energy  by  less  than  5%.  However,  the  calculated  interfaeial  energies  for 
stacking-faults  show  strong  dependence  on  the  adopted  supereell  size,  as  it  is  seen  from  Table  1 
(II  and  12)  that  the  calculated  interfaeial  energies  are  reduced  by  50%  when  the  supereell  size  is 
increased  from  8  to  32.  This  is  because  that  the  interfaeial  energies  for  the  II  and  12  stacking- 
faults  are  small  number  whose  calculations  have  reached  the  accuracy  limit  of  the  current  first- 
principles  method. 


We  assume  that  the  larger  the  supereell  sizes  are  the  more  aeeurate  of  the  ealeulated 
interfaeial  energies  are.  For  the  six  largest  supereells  of  the  six  types  of  defeets,  our  ealeulated 
interfaeial  energies  are  85.5,  81.0,  118.1,  120.0,  8.1,  and  21.8  mJ/m^,  respeetively,  for  the  (1011) 
mirror  refleetion,  the  (1011)  mirror  glide,  the  (1012)  mirror  refleetion,  the  (1012)  mirror  glide, 
the  II  staeking  fault,  and  12  stacking  fault.  In  overall,  these  numbers  are  in  good  agreements  with 
the  reported  results  from  the  previous  publications.  Yoo  et  al.  [14]  calculated  the  (1011)  and 
(1012)  mirror  reflection  twins  employing  the  ab  initio  method  (no  information  was  given  on  the 
employed  potential  and  the  approximation  to  the  exchange-correlation  energy  by  Yoo  et  ah),  and 
reported  the  calculated  interfaeial  energy  values  of  70  and  114  mJ/m  ,  respectively;  Morris  et  al. 
employed  the  same  ab  initio  method  (the  same  group  of  Yoo  et  al.  mentioned  above)  to  calculate 
the  (1012)  mirror  reflection  and  the  (1012)  mirror  glide  twins  using  20-atom  supereells  and 
reported  the  same  interfaeial  energy  of  114  mJ/m  for  both  the  (1012)  mirror  reflection  and  the 
(1012)  mirror  glide  twins.  The  calculated  interfaeial  energies  of  the  II  and  12  stacking  faults  by 
Smith  et  al.  [12]  were,  respectively,  9  and  18  mJ/m  using  the  24-atom  supereells  (12  layers  of 
Mg)  and  the  ABINIT  package  [19]  within  GGA.  The  calculated  interfaeial  energies  of  the  II  and 
12  stacking  faults  by  Che  tty  and  Weinert  [10]  were,  respectively,  10and22mJ/m  using  the  24- 
atom  supereells  (12  layers  of  Mg)  within  the  local  density  approximation  (LDA)  [20].  The 
calculated  interfaeial  energies  of  12  stacking  fault  by  Uesugi  et  al.  [13]  was  16  mJ/m  using 
CASTEP  package  [21]  within  GGA.  Other  results  for  stacking-fault  energies  are  from  the 
embedded  atom  method  (EAM)  by  Hu  et  al.  [4]  who  reported  the  interfaeial  energies  are  4,  and  8 
mJ/m  ,  respectively,  for  the  Hand  12  stacking  faults,  respectively. 

In  summary,  the  interfaeial  energies  for  four  types  of  twin  boundaries  and  two  types  of 
stacking  faults  observed  in  metal  magnesium  have  been  obtained  through  first-principles 


calculations.  It  therefore  offers  a  unified  pieture  of  the  interfaeial  energies  for  these  lattice 
mismatches  within  the  same  theoretical  framework.  The  caleulated  values  can  also  serve  as  the 
input  for  the  future  simulation  of  the  growth  process  of  these  planar  defeets  or  the  estimation  of 
the  effeets  of  these  lattice  mismatehes  on  the  mechanical  properties  of  magnesium  alloys. 
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Figure  1  (color  online).  The  (1011)  and  (1012)  planes  in  hep  metals.  The  heavy  (blue)  and 
light  (yellow)  balls  are  used  to  represent  the  A  and  B  atomic  layers  in  c  direction  conventionally 
used  for  the  hep  structure.  The  grey  box  shows  the  orthogonal  cell  mapped  from  the  hep  cell. 


Figure  2  (color  online).  Structures  of  (1012)  mirror  reflection  twin-boundary  (left  panel)  and 
(1012)  mirror  glide  twin-boundary  (right  panel).  See  also  Fig.  1  for  meanings  of  the  balls. 


Figure  3  (color  online).  Structures  of  (1011)  mirror  reflection  twin-boundary  (left  panel)  and 
(1011)  mirror  glide  twin-boundary  (right  panel).  See  also  Fig.  1  for  meanings  of  the  balls. 


Table  1.  Computational  settings  and  interfaeial  energies.  For  the  II  and  12  staeking  faults,  the 
reported  stacking  fault  energies  in  the  literatures  [6-8,  10,  12,  13]  have  been  divided  by  a  factor 
of  two  since  II  and  12  contain  two  interfaces  [12], 
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